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1 METHOD OF OBTAINING 

2 PORE PRESSURE AND FLUID SATURATION CHANGES IN 

3 SUBTERRANEAN RESERVOIRS BY FORWARD MODELING 

4 

5 BACKGROUND 

6 The use of seismic data in the analysis and modeling of subterranean reservoirs 



7 containing hydrocarbons and other fluids is known. Typically, such data are gathered 

8 through the use of a source of seismic energy and one or more receivers respectively 

9 located on a ground or water surface over a subterranean region of interest. The source 

10 is used to produce a seismic pulse, burst, or similar energy which travels generally 

1 1 downward and away from the source, into the subterranean material of the region under 

12 examination. 

13 As the seismic pulse encounters a change in material properties, most notably at 

14 an interface between one type of subterranean material and another, some of the 

15 seismic pulse energy is reflected back toward the surface. The receiver or receivers 

16 detect this reflected pulse energy and record corresponding data, often with respect to 

17 other parameters of interest such as linear distance from the particular receiver to the 

18 source, time-of-flight (i.e., time between emission of source pulse and detected 

19 reflection), amplitude of the detected reflection, angle of incidence of the detected 

20 reflection relative to the ground (or water) surface plane or some other datum, etc. 

21 Thus, the presence of the interface can be detected through later analysis of the 

22 detected and recorded pulse reflection data. 



23 Generally, such pulse reflection and associated parameter data have been used 

24 to model, or estimate, the depths of these subterranean material interfaces and to 

25 present this information in the form of a cross-sectional elevation plot of the 

26 subterranean region of interest. However, such a plot often fails to provide other 

27 desirable information regarding the present physical state of a subterranean reservoir 

28 containing hydrocarbons or other fluids. 

29 Therefore, it is desirable to provide a method and apparatus for modeling various 

30 other subterranean physical parameters, and to present that model in the form of planar 

31 view representation (as well as 3D view presentation) of the subterranean region of 

32 interest. 
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1 SUMMARY 

2 One embodiment of the present invention provides for a method of modeling 

3 seismic data. The method includes deriving a time-lapse data set from a first seismic 

4 data set and a second seismic data set, and deriving a forward-modeled time-lapse data 

5 set including a plurality of values. The method further includes sorting the plurality of 

6 values into a plurality of bins corresponding to the forward-modeled time-lapse data set, 

7 selecting a plurality of optimal values from the plurality of bins, and then mapping the 

8 plurality of optimal values using the time lapse data set. The method also includes 

9 calibrating the plurality of optimal values. The method further includes plotting the 

10 plurality of calibrated optimal values. 

11 Another embodiment provides for a method of modeling seismic data 

12 corresponding to a subterranean reservoir containing hydrocarbons. The method 

13 includes calibrating a first seismic data set and a second seismic data set, and then 

14 subtracting the calibrated second seismic data set from the calibrated first seismic data 

15 set to derive a time-lapse data set. The method further includes deriving a forward- 

16 modeled time-lapse data set including a plurality of physical parametric values, sorting 

17 the plurality of physical parametric values into a plurality of bins corresponding to the 

18 forward-modeled time-lapse data set, and selecting a plurality of optimal physical 

19 parametric values from the plurality of bins of physical parametric values. The method 

20 also includes mapping the plurality of optimal physical parametric values to a 

21 corresponding plurality of subterranean locations using the time-lapse data set, and 

22 calibrating the plurality of optimal physical parametric values. The method also includes 

23 plotting the plurality of calibrated optimal physical parametric values as a visual 

24 representation of the subterranean reservoir containing hydrocarbons. 

25 Yet another embodiment provides for a computer which includes a processor and 

26 a computer-readable storage medium coupled in data communication with the 

27 processor. The computer-readable storage medium stores a first data set and a second 

28 data set and a plurality of rock physics relationships and a program code. The program 

29 code is configured to cause the processor to calibrate each of the first and second data 

30 sets, and then to subtract the calibrated second data set from the calibrated first data set 

31 to derive a time-lapse data set. The program code is further configured to cause the 

32 processor to calculate a forward-modeled time-lapse data set including a plurality of 

33 parametric values using selected ones of the plurality of rock physics relationships. The 

34 program code is still further configured to sort the plurality of parametric values into a 

35 plurality of bins corresponding to the forward-modeled time-lapse data set, and to select 
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1 a plurality of optimal parametric values from the plurality of parametric values sorted into 

2 the plurality of bins. The program code is further configured to cause the processor to 

3 map the plurality of optimal parametric values to a corresponding plurality of 

4 subterranean locations using the time-lapse data set, calibrate the plurality of optimal 

5 parametric values, and to plot the plurality of calibrated optimal parametric values to 

6 visually represent at least one spatially distributed physical characteristic of a 

7 subterranean reservoir of hydrocarbons. 

8 These and other aspects and embodiments will now be described in detail with 

9 reference to the accompanying drawings, wherein: 
10 

11 DESCRIPTION OF THE DRAWINGS 

12 The patent or application file contains at least one drawing executed in color. 

13 Copies of this patent or patent application publication with color drawing(s) will be 

14 provided by the Office upon request and payment of the necessary fee. 

15 Fig. 1 is a side elevation sectional view depicting a field seismology arrangement 

16 in accordance with an embodiment of the present invention. 

17 Fig. 2 is a flowchart depicting a method of modeling seismic data in accordance 

18 with another embodiment of the present invention. 

19 Fig. 3 is a mapping diagram referring to exemplary data plots 3A-3D in 

20 accordance with the present invention. 

21 Figs. 3A through 3D are data plots depicting exemplary P-wave and S-wave 

22 seismic data collected at times T1 and T2 in accordance with the present invention. 

23 Fig. 4 is a mapping diagram referring to exemplary data plots 4A-4B in 

24 accordance with the present invention. 

25 Figs. 4A and 4B are data plots depicting exemplary inverted time-lapse P-wave 

26 and S-wave seismic data in accordance with the present invention. 

27 Fig. 5 is a block diagram depicting a data cube in accordance with an 

28 embodiment of the present invention. 

29 Fig. 6 is a block diagram depicting the data cube of Fig. 5 and a corresponding 

30 data array in accordance with an embodiment of the present invention. 

31 Fig. 7 is a locus plot in accordance with an embodiment of the present invention. 

32 Fig. 8 is a block diagram depicting data mapping in accordance with an 

33 embodiment of the present invention. 

34 Fig. 9 is a mapping diagram referring to exemplary data plots 9A-9D in 

35 accordance with the present invention. 
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1 Figs. 9A through 9D are data plots depicting exemplary calibrated and actual 

2 (true) saturation and pore pressure values in accordance with the present invention. 

3 Fig. 10 is a block diagram depicting a data acquisition and processing system in 

4 accordance with yet another embodiment of the present invention. 
5 

6 DETAILED DESCRIPTION 

7 In representative embodiments, the present teachings provide methods and 

8 apparatus for acquiring and processing seismic data corresponding to a subterranean 

9 region of interest, typically containing hydrocarbons, and to plot a final processed data 

10 set as a graphic representation of time-lapse changes to various selected physical 

1 1 parameters of the subterranean region of interest. 

12 Turning now to Fig. 1, a side elevation sectional view depicts a field seismology 

13 arrangement 100 in accordance with an embodiment of the present invention. The 

14 arrangement 100 includes a source (i.e., emitter) of seismic energy 102. The source 

15 102 is located on a ground surface 104 over a subterranean region 106. The source 

16 102 can be defined by any suitable apparatus capable of producing a seismic (acoustic) 

17 pulse or vibrational energy which is directed generally into the subterranean region 106. 

18 The field seismology arrangement 100 further includes a pair of seismic detectors 



19 108 and 110, respectively (also known as "receivers", as for example geophones, or 

20 hydrophones when the surface 104 is underwater). Each of the seismic detectors 108 

21 and 110 rests on the ground surface 104, and is spaced apart from the source 102 by an 

22 offset distance 1 12 or 1 14, respectively. Each of the seismic detectors 1 08 and 110 can 

23 be defined by any detection device suitable for detecting and recording seismic energy 

24 pulses or reflections (described in detail hereafter) that arrive at the detectors 108 and 

25 110 after passing through the subterranean region 106. 

26 The subterranean region 106 includes three different material strata designated 

27 as 116, 118 and 120, respectively. Each of the material strata 116, 118 and 120 is 

28 respectively defined by a depth D1, D2 and D3. Furthermore, the region 106 includes 

29 an interface 122 between the material strata 116 and 118, an interface 124 between the 

30 material strata 118 and 120, and an interface 126 between the material strata 120 and 

31 an underlying region 121. It is assumed that each of the strata 116, 118 and 120 

32 includes a respective average material incompressibility, fluid content (or lack thereof), 

33 and other physical parameters that distinguish it from the other respective material 

34 strata. 
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1 Typical operation of the field seismology arrangement 100 is as follows: the 

2 source 102 produces a source seismic pulse P1 of known amplitude A1. The source 

3 pulse P1 is directed into the subterranean region 106 and proceeds initially through the 

4 material strata 116, striking the interface 122 at angle of incidence AN1 . A portion of the 

5 energy of pulse P1 is reflected from the interface 122 back toward the surface 104, as 

6 reflection pulse P2, at an angle of reflection AN1. The pulse P2 arrives at surface 104 

7 within detectable vicinity to the seismic detector 108, having amplitude A2 upon arrival. 

8 The seismic detector 108 records data corresponding to the detection of the 

9 reflection pulse P2. This recorded data can include, for example, a detected amplitude 

10 corresponding to amplitude A2 of the pulse P2, the detected angle of reflection (i.e., 

11 incidence) AN1 of the pulse P2, the offset distance 112, the arrival time of the detected 

12 pulse P2 relative to the (known) time of emission of the source pulse P1, etc. It will be 

13 appreciated that the seismic detector 108 only transitorily "records" data, and that in fact 

14 the detector 108 transmits the data to a permanent recording station (not shown) for 

15 recording on computer readable media such as a magnetic tape or a hard disk drive. 

16 As the source pulse P1 continues into the subterranean region 106, similar 

17 reflection pulses P3 and P4 are reflected from interfaces 124 and 126, respectively. The 

18 reflection pulse P3 is initially reflected from the interface 124 back toward the surface 

19 104 at an angle of AN2, and is then refracted at the interface 122 to a new angle of 

20 incidence AN1. Thus, the reflection pulses P3 and P4 arrive at the surface 104 within 

21 detectable vicinity of the seismic detectors 108 and 110, respectively. The seismic 

22 detectors 1 08 and 1 1 0 then record data corresponding to the reflected pulses P3 and 

23 P4. This recorded data can include any or all of the various characteristics described 

24 above in regard to the pulse P2. 



25 The data thus received by the seismic detectors 108 and 110 are recorded and 

26 then communicated to a suitable analytical apparatus (i.e., a computer) for analysis by 

27 way of the method of the present invention, described in detail hereafter. 

28 It is to be understood that the field seismology arrangement 100 of Fig. 1 is 

29 intended to convey conceptual information regarding the acquisition of seismic data as 

30 used in the methods of the present invention, and that a similar field seismology 

31 arrangement (not shown) can include any suitable number of seismic energy sources 

32 (i.e., source 102) and seismic detectors (i.e., detector 108) arranged in a linear or matrix 

33 pattern on a ground surface (i.e., surface 104) or on a water surface in the case where 

34 the ground surface 104 is submerged under water. Thus, the field seismology 

35 arrangement 100 of Fig. 1 is exemplary of the data acquisition of the present invention 



Case FO01-P03 



1 and does not represent the only type of such arrangement that can be used, nor does it 

2 represent the only type of field conditions under which data acquisition can be performed 

3 within the context of the present invention. 

4 Fig. 2 is a flowchart depicting a method 200 in accordance with the present 

5 invention. Simultaneous reference shall be made, as directed, to Figs. 1 and 3 through 

6 9D, during the description of the method 200 of Fig. 2. 

7 In step 202, amplitude-versus-offset (hereafter, AVO) field seismic data are 

8 acquired through the use of field seismology (e.g., the field seismology arrangement 100 

9 of Fig. 1). These AVO data are acquired at a first time T1 and at a later time T2, and are 

10 generally collected over a subterranean region of interest such that an area at a given 

11 depth (i.e., stratum) is represented. AVO data for both generally planar and volumetric 

12 regions of interest can be suitably acquired. 

13 in step 204, the AVO data are inverted to seismic impedance data using standard 



14 mathematical techniques. As a typical AVO data set is relatively vast, such inversion is 

15 generally done by way of electronic computer (see Fig. 10). The inverted data sets 

16 include P-wave pseudo impedance data (hereafter, IP') and S-wave pseudo impedance 

17 data (hereafter, IS'), and are referred to as IP'(T1), IS'(T1), IP'(T2), and IS'(T2). 

18 In step 205, the P-wave and S-wave pseudo impedance data IP'(T1), IS'(T1), 

19 IP'(T2), and IS'(T2) are calibrated so as to correspond more closely with actual field 

20 conditions. This calibration can be performed in a number of different ways; non-limiting 

21 examples include: calibrating the pseudo values against like kinds of data (i.e., P-wave 

22 and S-wave impedance data) measured at selected well bores; or modeling the pseudo 

23 values using rock physics relationships (properties). Combinations of these or other 



24 calibration methods can also be used. 

25 It is to be understood that such calibration is not necessarily linear in nature. In 

26 any case, the calibration method yields calibrated P-wave and S-wave impedance data 

27 IP (T1), IS (T1), IP(T2), and IS(T2) for the subterranean region under consideration. 

28 Exemplary plots of such calibrated data IP(T1), IS(T1), IP(T2), and IS(T2) are 

29 respectively depicted in Figs. 3A, 3B, 3C, and 3D. Within Figs. 3A-3D, color is indicative 

30 of the signal or seismic wave impedance, with red indicative of seismic wave impedance 

31 values at the lower end of the represented scale, and blue used to indicate seismic 

32 impedance values at the higher end of the represented scale. The color white (or 

33 colorless) is used to indicate mid-scale seismic impedance values. 

34 As such, Figs. 3C and 3D each depict the presence of two spot locations where a 

35 generally distinct decrease in seismic wave impedance has occurred, relative to the 
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1 same locations in Figs. 3A and 3B. These spot decreases (i.e., changes) in seismic 

2 wave impedance correspond to physical changes within the subterranean region of 

3 interest, such as, for example, changes in porosity, pore pressure, saturation, etc. 

4 In step 206, the inverted AVO data matrices derived in step 204 are subtracted in 

5 accordance with the two following formulas, thus providing the indicated time-lapse data 

6 (i.e., recorded field data): 

7 1) TL(IP) = [IP(T2)-IP(T1)] 

8 2) TL(IS) = [IS(T2)-IS(T1)] 



9 It is important to note that time-lapse data required by this embodiment of the present 

10 invention need not be seismic impedance (i.e., inverted AVO) data as described above 

11 for step 204. Any time lapse data that can be forward-modeled from a set of physical 

12 parameters within the context of the present invention can potentially be used in step 

13 206. Thus, in another embodiment (not shown) of the present invention, a method using 

14 appropriate AVO data acquired in step 202 can dispense with (i.e., skip over) the 

15 inverting of step 204 above and proceed directly to the calibrating of step 205 above. 

16 Other methods (not shown) in accordance with other embodiments of the present 

17 invention can also be used. 

18 Exemplary plots of such time-lapse data are respectively depicted in Figs. 4A and 

19 4B. As described above, color is indicative of (i.e., in correspondence to) seismic wave 

20 impedance value within Figs. 4A and 4B. The time-lapse data TL(IP) and TL(IS) derived 

21 within step 206 result in a generally more pronounced and distinct indication of physical 

22 changes within the subterranean region. As shown in Figs. 4A and 4B, the two spot 

23 locations described above are clearly pronounced, and are assumed to be indicative of a 

24 subterranean physical change of interest (i.e., porosity, saturation, etc.), or a 



25 combination of such changes. 

26 In step 208, selected known rock physics relationships and corresponding 

27 formulas are used to compute forward-modeled time-lapse data (i.e., synthetic data) 

28 FMTL(lp) and FMTL(ls). Typically, these relationships include such physical parameters 

29 as pore pressure, fluid saturation, and rock porosity. Such calculations can be generally 

30 represented by the two following formulas: 

31 3) FMTL(IP) = F1 [TL(Saturation), TL(Pore Pressure), TL(Porosity)] 

32 4) FMTL(IS) = F2[TL(Saturation), TL(Pore Pressure), TL(Porosity)] 

33 where F1 and F2 are selected rock physics relationships (i.e., formulas) that are 

34 functions of the desired physical parameters of porosity, saturation, and pore pressure. 

35 Other physical parameters, by way of their associated formulas, can also be considered 
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1 such as, for example, temperature, salinity, gas-to-oil ratio (GOR), gas gravity, 

2 overburden pressure, etc. 

3 Reference is now made to Fig. 5. These rock physics calculations are generally 

4 used to construct a data cube 300 of the three exemplary parameter types, respectively 

5 depicted as TL pore pressure 302, TL saturation 304, and TL porosity 306. The data 

6 cube 300 thus includes a plurality of three-dimensional cells 308, which respectively 

7 contain the forward-modeled time-lapse data FMTL(IP) and FMTL(IS) pair values 

8 corresponding to the coordinates (i.e., parametric values) of the particular cell 308. An 

9 exemplary cell 310 of the plurality of cells 308 is depicted, which includes corresponding 
10 FMTL(IP) and FMTL(IS) data pair contents 312. The data cube 300 therefore 



1 1 represents a three-dimensional data model of the subterranean region of interest. It is to 

12 be understood that if other physical parameters are considered, a corresponding data 

13 cube (not shown) can have four or more dimensions. 

14 it is important to note that any given forward-modeled time-lapse data pair 

15 FMTL(IP) and FMTL(IS) can result from more than one corresponding set of physical 

16 parameters - that is, more than one cell 308 within the data cube 300. Typically, any 

17 given time-lapse data pair (for example, data pair 312) results from several 

18 corresponding sets of physical parameters, which can be visualized as rays or arcs of 

19 adjacent or near-adjacent, associated cells 308 within the data cube 300. 

20 In step 210, the forward-modeled time-lapse physics data within data cube 300 is 



21 sorted. Reference is now made to Fig. 6. A two-dimensional array 314 including a 

22 plurality of data bins 316 is constructed, with each data bin 316 defined by coordinates 

23 corresponding to the values of a particular forward-modeled time-lapse data pair (i.e., 

24 FMTL(IP) and FMTL(IS)). As described above, any particular forward-modeled time- 

25 lapse data pair can correspond to a plurality of derived physical parameters, and 

26 therefore several associated physical parameter sets, or vectors, can be sorted into any 

27 particular data bin 316 of the array 314. As depicted in Fig. 6, an exemplary data bin 

28 318 of the plurality of data bins 316 is depicted, which contains (i.e., includes) an 

29 associated physical parameter vector 320. As depicted in Fig. 6, Pp' refers to pseudo 

30 pore pressure, Sw' refers to pseudo water saturation, and cp (phi) refers to porosity. 



31 The sorting process is conducted in an exhaustive fashion until all the physical 

32 parameter vectors (i.e., 320) have been sorted into their respective data bins 316 within 

33 the array 314. 

34 In step 212, the contents of each data bin 316 within the array 314 are compared 

35 (i.e., searched) to a predetermined, selected parameter value, so as to determine which 
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1 particular physical parameter vector represents the "optimal" such vector within each 

2 data bin 316. For example, one approach for conducting this search is to compare each 

3 of the physical parameter vectors with an average or sample porosity value for the 

4 subterranean region under consideration. This comparison value can be predetermined, 

5 say, by use of appropriate field instrumentation deployed within a borehole or similar 

6 arrangement (not shown). Other search and comparison techniques can be used. 

7 Reference is now made to Fig. 7, which is a plot 330 depicting the contents of the 

8 data bins 316 of the array 314, and is provided to assist in an understanding of the 

9 optimal value search operation of step 212 of the method 200. The plot 330 is formatted 

10 with a horizontal axis scaled to represent time-lapse pressure values 332, and a vertical 

11 axis scaled to represent time-lapse saturation values 334. Each of the data bins 316 of 

12 the array 314 has its physical parameter vectors plotted as a single locus 336 of values 

13 on the plot 330 (that is, there is one plotted locus 336 for each data bin 316). 

14 Within each locus 336 is a selected optimum parameter pair (i.e., vector) value 

15 338, including corresponding pressure 332 and saturation 334 values, as determined by 

16 the comparative search described above. The optimum pressure 332 and saturation 

17 334 parameter pairs 338 are extracted for further use as described hereafter. 

18 In step 214, the optimum parameter pairs 338 are mapped to their corresponding 

19 locations within the subterranean area under consideration. Steps 210-214 are 

20 generally referred to as inversion. Reference is now made to Fig. 8, which depicts a 

21 parameter mapping schema 350 in accordance with the present invention. To begin, the 

22 time-lapse data pair TL(IP) and TL(IS) from step 206 for each location within the 

23 subterranean region under consideration is isolated, one data pair at a time. The data 

24 bin 316 within the array 314 that corresponds to the values of an isolated time-lapse 

25 data pair is then referenced, and its parameter vector contents considered. The optimal 

26 parameter vector 338 within that data bin 316, as determined in step 212 above, is then 

27 extracted from the data bin 31 6. 

28 The desired discrete physical parameters within the optimal parameter vector 338 

29 are then associated with the subterranean location of the original time-lapse data pair 

30 TL(IP) and TL(IS). As depicted in Fig. 8, the time-lapse data pair 312 is associated with 

31 a location 352 within the subterranean region corresponding to the original AVO data. 

32 Therefore, the particular optimal parameter vector 338, depicted as a vector 354, is also 

33 associated with the same location 352. As particularly depicted in Fig. 8, a pore 

34 pressure parameter 356 and a saturation parameter 358 of the vector 354 are 

35 associated with the location 352. These mapped, optimal parameters (i.e., pore 
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1 pressure 356 and saturation 358) are also referred to as pseudo values, and are 

2 designated in Fig. 8 as TL Press* and TL Sat', respectively. 

3 The mapping process of step 214 is generally repeated as described above, until 

4 optimal physical parameters are associated with each location within the subterranean 

5 region corresponding to the original AVO data. 

6 In step 216, the pseudo values (i.e., TL Press 1 356 and TL Sat' 358) mapped in 

7 step 214 above are calibrated so as to correspond more closely with actual field 

8 conditions. This calibration can be performed in a number of different ways; non-limiting 

9 examples include: calibrating the pseudo values against like kinds of data (i.e., pore 

10 pressures and saturations) measured at selected well bores; calibrating the pseudo 

1 1 values against a flow model of the subterranean region of consideration; or modeling the 

12 pseudo values against rock physics relationships (properties) in which only pore 

13 pressure changes or saturation changes. Combinations of these or other calibration 

14 methods can also be used. It is to be understood that such calibration is not necessarily 

15 linear in nature. In any case, the calibration method yields calibrated pore pressure and 

16 saturation data for the subterranean region under consideration. 

17 In step 218, the calibrated data from step 216 above are plotted to provide a 

18 2-dimensional representation of the subterranean region under consideration. 

19 Reference is now made to Figs. 9A-9D. This plot represents the time-lapse change in 

20 the physical parameters derived and calibrated as described above in regard to steps 

21 202-216 of the method 200. Figs. 9A and 9B represent the calibrated saturation and 

22 pore pressure data from step 216 above, respectively. Figs. 9C and 9D represent the 

23 actual (true) saturation and pore pressure data, respectively. Once again, the colors red 

24 and blue are used to indicate parameter value within the Figs. 9A-9D. Once the data 

25 plotting of step 218 is performed, performance of the method 200 is complete. 

26 It is to be understood that the method 200 of Fig. 2 represents one embodiment 

27 of the present invention, and that other methods (not shown) corresponding to other 

28 embodiments of the present invention can also be used. For example, the methods and 

29 teachings of the present invention can be used with other kinds of AVO data that can be 

30 forward-modeled from a set of physical parameters. The impedance data exemplified in 

31 method 200 represents just one of several possible approaches. As described above, 

32 under another embodiment of the present invention, for example, the inversion of step 

33 204 described above would be optional. 

34 Furthermore, other embodiments of the present invention can provide 

35 corresponding methods in which the certain steps or operations are performed 
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1 substantially in parallel with (i.e., concurrent to) other certain steps. For example, 

2 another embodiment (not shown; see Fig. 2) can provide for the performing of steps 202 

3 through 206 above substantially in parallel with the performing of steps 208 through 212 

4 above, wherein the respective results of these substantially parallel operations (e.g., the 

5 time-lapse data set TL(IP) and TL(IS), and the optimum parameter pairs 338) are then 

6 used to perform steps 214 through 218. Other methods and other embodiments of the 

7 present invention are also possible. 

8 Fig. 10 is a block diagram depicting a data acquisition and processing system 

9 (hereafter, system) 400 in accordance with yet another embodiment of the present 

10 invention. The system 400 includes a field seismology arrangement 402. As depicted in 

11 Fig. 10, the field seismology arrangement 402 includes a plurality of seismic detectors 

12 (i.e., receivers) 408, 410, and 412, respectively. The field seismology arrangement 402 

13 is assumed to include at least one source of seismic energy (not shown), and any other 

14 elements as required and configured to acquire desired seismic (i.e., AVO or AVA, 

15 where AVA is amplitude-versus-angle type seismic data) data corresponding to a 

16 subterranean region 414 underlying the field seismology arrangement 402. The field 

17 seismology arrangement 402 is further configured to provide one or more acquired 

18 seismic data bundles 416 for processing with the balance of the system 400 as 

19 described hereafter. 

20 The system 400 also includes a computer 418. The computer 418 includes a 

21 processor 420 coupled in data communication with a computer-accessible memory 422. 

22 The memory 422 stores a first seismic data set 424 and a second seismic data set 426. 

23 The first seismic data set 424 is assumed to be received by the computer 418 and 

24 stored in the memory 422, prior to the computer 418 receiving and storing the second 

25 seismic data set 426. It will be appreciated that the data sets 424 and 426 can also be 

26 stored in a remote memory device which is accessible by the computer 418. 

27 Both the first and second seismic data sets 424 and 426 are delivered to the 

28 computer 418 as corresponding seismic data bundles 416, and can be delivered to the 

29 computer 41 8 by way of any satisfactory means. Non-limiting examples of such delivery 

30 means (not shown) can include data cable coupling, transferal by way of optical or 

31 magnetic storage media, radio telemetry linking, etc. Those of skill in the 

32 instrumentation and related arts can appreciate that any number of satisfactory seismic 

33 data 416 delivery means can be utilized within the scope of the present invention, and 

34 that further elaboration is not required for purposes herein. 
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1 The memory 422 further stores a program code 428 that is executable by the 

2 processor 420. The program code 428 is configured to cause to the processor 420 to 

3 substantially perform the method 200 of Fig. 2 as described above. The memory 422 

4 also stores a plurality of rock physics relationships 430, which are selectively accessed 

5 and used by the processor 420 during execution of the program code 428. 

6 The system 400 also includes a monitor 432 that is coupled in signal 

7 communication with the computer 418. The monitor is configured to provide a user 

8 visible data plot 434 under the control of the processor 420 during execution of the 

9 program code 428. The system 400 further includes a printer 436 coupled in signal 

10 communication with the computer 418. The printer 418 is configured to provide a 

11 hardcopy data plot 438 under the control of the processor 420 during execution of the 

12 program code 428. 

13 The computer 41 8 is further understood to include a plurality of other elements as 



14 desired and/or required for normal operation, which are not shown in Fig. 10. Such 

15 elements (not shown) can include, for example, a user keyboard, a user mouse, a power 

16 supply, etc. One of skill in the computing arts can appreciate that such elements can be 

17 respectively included with the computer 418 and configured as desired, and that further 

18 elaboration is not required for an understanding of the present invention. 

19 Typical normal operation of the system 400 is as follows: The field seismology 

20 arrangement 402 acquires the first seismic data set 424, and at some predetermined 

21 period of time thereafter, the field seismology arrangement 402 acquires the second 

22 seismic data set 426. The first and second seismic data sets 424 and 426 are delivered 

23 to the computer 418 as respective seismic data bundles 416, which stores them 

24 accordingly within the memory 422. 

25 Next, execution of the program code 428 by the processor 420 is initiated by a 

26 user (for example, by way of a user keyboard or mouse, not shown). The program code 

27 428 then causes the processor 420 to selectively access the first and second seismic 

28 data sets 424 and 426, as well as the rock physics relationships 430, which are 

29 respectively stored in the memory 422. The processor 420 then uses the data sets 424 

30 and 426 and the rock physics relationships 430 to carry out (i.e., perform) the method 

31 200 of Fig. 2 substantially as described above, thus deriving a calibrated physical 

32 parameter data set corresponding to the first and second seismic data sets 424 and 426 

33 provided by the field seismology arrangement 402. 

34 The program code 428 then causes the processor 420 to plot the calibrated 

35 physical parameter data set using the monitor 432 and/or the printer 436, resulting in the 
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1 visible data plot 434 and/or the hardcopy data plot 438, respectively. The plot 434 

2 and/or 438 thus provides a visible representation of the selected time-lapse physical 

3 characteristics (i.e., porosity, pressure and/or saturation, etc.) of the subterranean region 

4 414. 

5 In this way, the system 400 of Fig. 10 provides a substantially automated data 

6 acquisition and processing system that performs the method of the present invention 

7 and provides resulting 2-dimensional data plots 434 and/or 438. The system 400 is 

8 particularly suitable for use in monitoring and analyzing time-lapse changes in various 

9 physical parameters of subterranean regions (i.e., region 414) that contain hydrocarbons 

10 such as crude oil, natural gas, etc, or other fluids such as water. It is to be understood 

11 that three-dimensional volumetric plots (not shown) corresponding to a case of three- 

12 dimensional inversion can also be provided under the present invention. 

13 Furthermore, it is to be understood that while the methods of the present 

14 invention described above consider first and second AVO data sets, any number of 

15 suitable data sets can also be considered within corresponding other embodiments (not 

16 shown) of the present invention. Within such embodiments (not shown), the methods 

17 and teachings of the present invention would typically be applied to any two suitable 

18 data sets at a time. 

19 While the above methods and apparatus have been described in language more 

20 or less specific as to structural and methodical features, it is to be understood, however, 

21 that they are not limited to the specific features shown and described, since the means 

22 herein disclosed comprise preferred forms of putting the invention into effect. The 

23 methods and apparatus are, therefore, claimed in any of their forms or modifications 

24 within the proper scope of the appended claims appropriately interpreted in accordance 

25 with the doctrine of equivalents. 
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